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A simple picture describes the results of recent treatments of partially- 
condensed, dilute, trapped Bose gases at temperature T > 0. The condensate 
wavefunction is nearly identical to that of a T = condensate with the 
same number of condensate atoms, Nq. The cloud of non-condensed atoms 
is described by the statistical mechanics of an ideal Bose gas in the com- 
bined potentials of the magnetic trap and the cloud-condensate interaction. 
We provide a physical motivation for this result, show how it emerges in 
the Hartree-Fock-Bogoliubov-Popov approximation, and explore some of its 
implications for future experiments. 

PA CS numbers: 3. 75. Fi, 67. 4 0. Db, 67. 90. +z 

1. Introduction 

Most recent experimentsBUi on dilute, magnetically-trapped, alkali- 
atom Bose gases have viewed phenomena which are well described by the 
zero-temperature (T = 0) mean-field theory of the Bose-Einstein conden- 
sate (BEC), in which virtually all the gas in the system resides in the con- 
densed state. The technique of evaporative coolingja used in all such ex- 
periments, grows the BEC by selective extraction of the non-condensed, 
"thermal," component of the gas, which is located at the outer edges of the 
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trap la The T = mean-field theory has been found to give a good account 
of many BEC properties observed in these systems. B0 New experiments! 1 ! 
have started to probe BECs at T > 0, and so the testing of alternative 
finite-temperature BEC theories has begun. 

This paper draws attention to common features emerging from sev- 
eral independent finite-temperature theories, which suggest that a relatively 
simple picture, which we call the "two-gas model," describes many of the 
properties of a dilute Bose gas with repulsive atomic pair interactions (scat- 
tering length a > 0). The two gases concerned are the condensate gas, the 
intrinsic properties of which are essentially independent of temperature; and 
the thermal, non-condensed gas, which behaves much like an ideal Bose gas 
at temperature T in an effective potential created by the condensate. This 
model emerges naturallv aaa limiting case of the Hartree-Fock-Bogoliubov- 
Popov (HFB-Popov)li3'EJ't3 and Hartree-Fock (HF)EL3 approximations, but 
its features seem also to be manifest in recent quantum Monte CarloE-3 and 
semiclassicalEl calculations. The two-gas picture offers some straightforward 
implications for interpretation of experiments and for further development 
of first-principles theories. 

Our evidence for the validity of this picture first emerged from large- 
scale numerical calculations, but its origin can be traced back qualitatively 
within the structure of finite-temperature field theory. In Sec. 2, we show 
how such a theory can plausibly lead to a two-gas scenario. Section 3 ex- 
plores some of the implications of the model. 



2. Two gas model as a limit of the HFB-Popov approximation 

The HFB-Popov equations have been derived elsewhere^ and we merely 
state the basic equations here. In the Heisenberg equation of motion, the 
Bose field operator, ^(r) is decomposed into a c-number wavefunction ipo(r) 
that describes a condensate of iVo atoms, and an operator ip(r) describing 
the non-condensate atoms: VK r ) = v / ^oV ; o( r ) + "0(r). In the HFB-Popov 
approximation, the wave function for a condensate of trapped atoms satisfies 
a generalized Gross-Pitaevskii (GP) equation: 



{h + Uq [No \Mr)\ 2 + 2n(r)] } ^(r) = M)(r), 



(1) 



where Hq = w V 2 + Vt rap (r) is the Hamiltonian for a single atom of mass M 
and position coordinate r; the trapping potential (for cylindrically symmet- 
ric systems of current interest) is given by T4rap(r) = M (u 2 p p 2 + ^ 2 z 2 ^j /2, 
with Up and co z the radial and axial angular frequencies of the trap; 

expresses the binary interaction between atoms; the chem- 
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ical potential /i, interpreted as the work required to add one more atom to 
the condensate, is treated as an eigenvalue; and t/^o(r) is normalized to unity. 
The function ra(r) is the density of the non-condensed component of the 

gas, 

n(r) = £ { [Mr)| 2 + \ Vj (r)\ 2 } N, + \ Vj (r)\ 2 } , (2) 



where 



Nj = [exp(PE j )-l}- 1 (3) 



is the Bose-Einstein factor, = k^T and fee is the Boltzmann constant. 
The quasi-particle excitation energies Ej and amplitudes Uj(r),Vj(r) are 
obtained by solution of the coupled HFB-Popov equations: 

£ Uj (r) + N U |V>o(r)| 2 ^(r) = E jUj (r), (4) 

Cv 3 {y) + N U \Mr)\ 2 Uj(r) = -E jVj (r), (5) 

where C = H$ + 2Uo<n(r) — (j,, and n(r) = iVo | C r ) 1 2 + ^( r ) i s the total 
trapped-atom density. 

In simple physical terms, Eqs. describe a condensate subject to 

interaction with itself and a thermal cloud, with the cloud being generated 
by thermal excitations of condensate quasi-particles (There is also a non- 
thermal contribution to this cloud, the so-called "quantum depletion" term 
represented by the rightmost term of Eq. (^) , but it is much smaller than the 
thermal component except near T = 0). To solve these equations for a given 
atomic species and trap configuration, we fix the values of T and A^o, and 
then determine all other quantities self-consistently, eventually obtaining 
the total number of trapped atoms, N, via 



N 



J dv n(r) = N + J2 N J (6) 



By carrying out a sequence of such (laborious) calculations, we can map 
out the {N.,Nq,T} phase diagram of the interacting Bose gas. We present 
elsewhere£3 a detailed comparison of the results of this approach with ex- 
perimental data for the 87 Rb condensate at JILA;! for the temperature- 
dependent quasiparticle excitation energies, HFB-Popov agrees with experi- 
ment to within 5% for temperatures from zero up to 65% of the temperature 
To of the phase transition for the corresponding ideal gas (corresponding to 
thermal fractions from less than 1% to about 50%). Although at present 
there are considerable discrepancies as T — > To, it seems that HFB-Popov is 
a useful working theory over a significant temperature range. 
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Several calculationallHl3 have shown that, for current experiments, the 
quantum depletion of a small condensate is negligible, i. e., J dr n(r)|y = o = 
J2j I dr \ vj(r)\ 2 « No. This justifies use of the approximation 

Vj(r) = 0, (7) 

which i s e quivalent to the Hartree-Fock approximation used by other au- 
thorsJl3EJ If we apply this approximation to Eq. (Q) (and neglect the con- 
tribution of n(r) to n(r)), we obtain an ordinary Schrodinger equation for 
Uj(r), 

h 2 



2M* + V «® 



E jUj (r), (8) 



where the effective potential, given by 

K ff (r) = y trap (r) + 2N U \Mr)\ 2 , (9) 

is that of the trap modified by the repulsive pair interaction between the 
thermally-excited atoms and the condensate density. 

If we consider the case of a relatively small thermal fraction, then we 
expect to find the condensate localized near the center of the trap, so that 
Eq. (H) presents the thermal cloud with a trap and repulsive core. Thus, at 
least the low-energy quasi-particle amplitudes Uj(r) will be expelled from 
the core, i. e., they will have little overlap with the condensate wavefunction. 
This then gives consistency of Eq. (||) with our initial approximation, Eq. 
(0). It also justifies the approximation that completes our portrayal of a 
two-gas system: because of the expulsion of quasi-particle amplitudes from 
the condensate, we assume that n(r) can be neglected in Eq.(|l]), so that the 
condensate wavefunction is determined by solving 



|^V 2 + y trap (r) + AWo|V<o(r)| 2 



V>o(r) = M)(r), (10) 



which is just the usual GP equation for a condensate of Nq atoms at T = 0. 

Thus, these arguments have led us to a simple picture in which the 
finite-temperature Bose system appears to be composed of two distinct 
gases. One of these gases, the BEC, is always effectively at zero tempera- 
ture, and is described by an equation which depends only on its own atomic 
population, No, and the trap parameters. The other gas, the thermal cloud, 
behaves as a normal Bose gas at finite temperature, sensing the presence 
of the condensate through an elastic interaction; it does not undergo Bose- 
Einstein condensation itself, but serves as an atomic reservoir for the BEC. 
This resembles the phenomenon of BEC of an ideal gas in an external po- 
tential, except that we account for interactions in the identification of the 
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ground state, and in the modification by the condensate of the external 
potential exposed to the thermal cloud. 



3. Implications of the two-gas model 

The two-gas model provides us with a straightforward way of computing 
the phase diagram of the dilute Bose gas for T <Tq. If there are A^o atoms 
in the BEC, we solve Eq. (|To| ) to obtain what we will call an equivalent zero- 
temperature condensate (EZC), i. e. the corresponding T = condensate 
that contains A"o atoms. The EZC solution provides us an orbital tpo(r;No) 
and chemical potential fi(No); with these in hand, we can construct V e s (r) 
and find the spectrum of Eq.(||). This procedure, which is independent of 
T, gives us the information we need to compute the equilibrium value of N 
for given values of Ao and T: we evaluate Eq. (||) from Eq. (|3|). 

In short, the EZC provides a "reference condensate" which, for a given 
set of trap parameters, describes all systems with the same number of con- 
densate atoms Nq. As we have shown elsewhere,t3 this model provides good 
agreement with the results of full HFB-Popov calculations of condensate and 
thermal densities and the critical temperature To; the emergence of an EZC 
can also be seen in the analysis by KrauthB of the results of his quantum 
Monte Carlo calculations, and in the recent quasi-classical calculations of 
Minguzzi et altB Comparison of the EZC condensate densities with those of 
HFB-Popov calculations is made in Fig. (jl|); this shows that even in cases 
where the condensate fraction / = Nq/N is as small as 0.1, the condensate 
is relatively unperturbed by the presence of the thermal cloud. 

Another straightforward implication of the two-gas model concerns the 
density profile of the thermal cloud. If we entertain the simple hypothesis 
that the cloud would be described by classical statistical mechanics of an 
ideal gas, then its density n(r) would be proportional to exp [— (3V e fi(r)]. 
Since V e a(r) is repulsive at small |r| and confining at large |r|, n(r) will 
attain its maximum awayfrom the center of the condensate. If we consider 
the Thomas-Fermi limitE-0 appropriate to large condensates, then for the 
case of a spherical condensate [14 rap (r) = Mu 2 r 2 /2] we find that 



V eS (r) 



Mf(2rl F -r 2 ), r < r TF 
^r 2 , ' r>r TF 



1/5 



where the Thomas-Fermi radius, ttf = 4^m^ ' > defines the sharp 
boundary of the condensate. Thus in this limit, n(r) is largest at the sur- 
face of the condensate, and its distribution becomes more localized as Nq 
increases, albeit slowly. 
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Fig. 1. Condensate densities in the z = plane for systems of Rb atoms in 
the JILA TOP trapta with radial frequency v p = 74 Hz. The solid line shows 
the condensate density as computed in the HFB-Popov approximation, and 
the dashed line is the density of the corresponding EZC. Case (a): A system 
of N = 13150 atoms at T = 70 nK, corresponding to / = 0.54 in the HFB- 
Popov approximation, i. e. Nq = 7106. (b) A system of N = 2000 atoms 
at T = 51 nK, corresponding to / = 0.1. The EZC densities are seen to be 
very close to those of the HFB-Popov approximation. 

The key qualitative aspects of this classical description are applicable 
to the quantum system, as shown in Fig.(^): this displays results of a full 
quantum-mechanical finite-number description, without any of the semiclas- 
sical continuous spectrum approximations made by other authors.li3'E3 This 
figure clearly suggests that quantitative interpretations of experimental data 
on finite-temperature condensates, e.g. determination of a condensate frac- 
tion from density measurements, will have to invoke some detailed model 
of the thermal distribution, since this distribution is neither monotonic nor 
close to the results obtained for a noninteracting gas. On the other hand, our 
current model suggests that the condensate and thermal densities are rela- 
tively distinctly segregated within the cloud, which may considerably sim- 
plify the qualitative understanding of some properties. Since it originates 
in the distinction between interactions of condensate and non-condensate 
atoms, this spatial segregation of the two gases seems to be a fundamental 
aspect of the behavior of inhomogeneous Bose gases, such as the trapped- 
atom systems of current interest. This may have interesting consequences 
for applications: for example, it may be possible to selectively extract con- 
densate vs. thermal atoms from a trap by appropriate positioning of a probe, 
thus obtaining an outcoupled matter wave with higher coherence than would 
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be otherwise expected. El In homogeneous systems, on the other hand, con- 
densate and non-condensate populations are intertwinedj_this is one of the 
essential features of the two-fluid model of liquid helium.E3 

Taking this idea further, we suggest that any property of a finite- 
temperature BEC should be compared in the first instance to that of the cor- 
responding EZC. In the two-gas model, we expect most of the T-dependence 
of a given quantity to be reduced to Ao-dependence. For example, in Fig. 
(jH) we show the quasiparticle excitation frequencies for the JILA TOP trap, 
over a range of temperatures relevant to recent experiments, as computed 
in the full HFB-Popov approximation and in the two-gas model. It is seen 
that the two methods agree up to temperatures quite close to the phase 
transition, so the main effect of finite temperature is renormalization of the 
value of Nq. An analogous result was seen in earlier calculations!^ cl for 
the homogeneous Bose-condensed gas of the temperature dependence of the 
speed of sound, which found it to be given by an equivalent T = expression 
adjusted for the temperature-dependence of the condensate density. 



4. Conclusions 

We find that condensate and thermal populations of a partially Bose- 
Einstein-condensed trapped-atom system separate out to a considerable ex- 
tent. Treating the condensate as uncoupled from the thermal cloud, and the 
thermal cloud as interacting with a static condensate potential, yields results 
similar to those that come from involved self-consistent field calculations. 
These results motivate the identification of the equivalent zero-temperature 
condensate (EZC) as a consolidating feature of finite-temperature systems. 
In this model, the main effect of finite temperature on the condensate is de- 
pletion of the condensate number. Condensate properties that depend only 
weakly upon No, such as the quasi-particle spectrum in the large-JVo limit 
(corresponding to the excitation frequencies of large condensates) ,0 should 
exhibit only weak temperature dependence. 



ACKNOWLEDGMENTS 

This work was supported in part by the U. S. National Science Foun- 
dation under grants PHY-9601261 and PHY-9612728, the U. S. Office of 
Naval Research, by MCS grant No. PAN/NIST-93-156, and by the U. K. 
Engineering and Physical Sciences Research Council. 



REFERENCES 



R. J. Dodd et al. 




2 4 6 8 10 

r/d 



Fig. 2. Thermal density for N ~ 2000 87 Rb atoms in a spherical trap with 
v r = uj/{2it) = 200 Hz, T = 75 nK, corresponding to / = 0.5 in the HFB- 
Popov approximation. The radius r is given in units of the characteristic 
length of the oscillator d ==yjfi/ {Moj). These parameters were previously 
used by Hutchinson et al.,EJ whose numerical HFB-Popov results we have 
reproduced and use here. The main figure shows the thermal density as 
computed in various approximations. Chain-dashed line: the confined ideal 
quantum gas; dashed line: full solution of the HFB-Popov equations; solid 
line: full quantum-statistical distribution as computed from two-gas model, 
which obtains N = 1965 (vs. the exact value of 2000) from iVo = 1000 
and T = 75 nK; dotted line: result of classical statistical mechanics applied 
to the two-gas model, using a fit to force N = 2000. The inset compares 
the HFB-Popov thermal density (solid line) with that of the condensate 
(dashed). Thus, even when the system is only 50% condensate, the peak 
condensate density is clearly much higher than that of the thermal cloud. 
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Fig. 3. HFB-Popov excitation frequencies (filled circles) for the m = 
(top), m = 2 (middle), and the m = 1 modes (bottom) for the JILA TOP 
trap with radial frequency v r = 129 Hz and N = 2000 87 Rb atoms, vs. 
temperature in units of To. Overlaid (solid lines) are the frequencies for a 
zero-temperature system with the same number Nq of condensate atoms as 
in the finite-temperature system. 
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